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Stimulated emission of photons from a large number of atoms into the mode of a strong 
■ Wght field is the principle mechanism for lasing in "classical" lasers 1 . The onset of lasing is 

CO 

marked by a threshold which can be characterised by a sharp increase in photon flux as a 

O . 

function of external pumping strength 2 3 . The same is not necessarily true for the fundamen- 



tal building block of a laser: a single trapped atom interacting with a single optical radiation 
models!. It has been shown that such a "quantum" laser can exhibit thresholdless lasing 
in the regime of strong coupling between atom and radiation field 7 8 . However, although 
theoretically predicted 9 12 , a threshold at the single-atom level could not be experimentally 
observed so far. Here, we demonstrate and characterise a single-atom laser with and without 
threshold behaviour by changing the strength of atom-light field coupling. We observe the 
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establishment of a laser threshold through the accumulation of photons in the optical mode 
even for a mean photon number substantially lower than for the classical case. Furthermore, 
self-quenching occurs for very strong external pumping and constitutes an intrinsic limita- 
tion of single-atom lasers Moreover, we find that the statistical properties of the emitted 
light can be adjusted for weak external pumping, from the quantum to the classical domain. 
Our observations mark an important step towards fundamental understanding of laser op- 
eration in the few-atom limit 13 including systems based on semiconductor quantum dots ^ 
or molecules^. 

Conventional lasers reach their threshold when the rate at which photons are coherently fed 
into the cavity mode exceeds losses and stimulated emission dominates over spontaneous emission 
into the lasing mode 2 3 . The most prominent characteristic of this threshold is a steep increase 
of photon flux from the laser while increasing the strength of the external pump exciting the gain 
medium. At the same time intensity fluctuations are reduced-^ 12 * 16 * 17 ! These characteristics are 
induced by the non-linear process of photon stimulated emission of photons into the lasing mode. 
At the single-atom level, characteristics that strongly deviate from standard laser properties have 
been predicte d 1 1 * 1 2 * 1 6 * 1 in particular regarding the question of the existence of a threshold. For 
strong atom-cavity coupling no threshold exists. In this regime, photons are coherently scattered 
into the cavity mode through stimulation by the vacuum field ®, More importantly, strong cou- 
pling introduces a photon-blockade, which suppresses cavity mode populations with more than 
one photorpjl. As a consequence, photon-stimulated emission enabling the threshold is absent, the 
laser threshold disappears and non-classical radiation is emitted^- 9 * 12 * 21 1 It was shown theoretically 
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that a one-atom laser can exhibit a threshold 11 ^ in an intermediate coupling regime. Here we ex- 
plore a single-ion device in this regime that exhibits a laser threshold for strong external pumping 
for which photon- stimulated processes start playing a role, owing to the non-zero population of 
cavity modes with more than one photon. 

The schematic experimental setup of our ion trap cavity-QED system is shown in Fig. la. 
Using frequency and intensity stabilised lasers, we excite a single 40 Ca + ion which is confined in 
a linear Paul trap at the centre of a near-concentric optical cavity^. The cavity provides a coherent 
coupling rate between the ion and a selected optical mode of g/2n = 1.3(1) MHz, limited by the 
localization of the ion along the cavity standing wave. With a cavity finesse of 70000 at 866 nm, 
the field decay rate from the cavity is given by kIItx = 54(1) kHz. A Hanbury-Brown&Twiss setup 
consisting of a beamsplitter and two avalanche photodiodes is used to analyse photons leaving the 
cavity. It provides us with the steady-state population of photons in the cavity 22 , n ss , and allows us 
to compute the second order correlation function, g^(r) (see supplementary information). 

In Fig. lb we present a reduced three level picture of the 40 Ca + ion's electronic structure 
which allows us to qualitatively describe the excitation process. Using an off-resonant drive laser, 
we achieve an effective coupling on a vacuum- stimulated Raman transition 4 between the electronic 
ground state S1/2 and the metastable state D 3 / 2 . A transfer from S1/2 to D 3 / 2 deposits a photon in 
the cavity. The effective coupling between these levels is g c s ~ g ■ si, where we have defined 
the drive parameter si = fii/2|Ai|, with denoting the Rabi frequency of the 397 nm driving 
laser and Ai its fixed detuning from the Si/ 2 - P1/2 resonance. To close the excitation cycle, an 
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infrared recycling laser at 866 nm, with a resonant Rabi frequency f2 2 and detuning A 2 , transfers 
population from the D 3 / 2 to the Pi/ 2 state. From there, spontaneous emission projects the ion back 
to the S1/2 ground state at a rate (27i/27r) = 20 MHz. The corresponding recycling rate is given 
by 27 r = 27! ,§2 with the recycling parameter s 2 = f2 2 /2A 2 . Here, we define the effective detuning 
A 2 = y/A 2 + (fi 2 /2) 2 + 7 2 where 27 = 2(71 + 72) is the total decay rate from the Pi/ 2 state 
and 27 2 /27r = 1.69 MHz is the decay rate back to the D 3 / 2 state. The total cycling rate r tot for 
external pumping, i.e. the rate at which one photon is added to the cavity mode, is thus given 
by l/r tot ~ l/2g cS + l/2j r . The strength of external pumping is then controlled by the drive 
and recycling lasers; however, note that the drive parameter s\ is kept constant in our experiments 
where the external pumping strength is only varied by the recycling parameter. Furthermore, the 
drive and recycling lasers introduce the broadenings 7^ ~ 7s 2 and 7^ ~ 7s 2 , of the Si/ 2 and 
D3/2 states, respectively. The total linewidth of the Si/ 2 <H- D 3 / 2 Raman transition is therefore 
given by 27 eff = 2(7^ + 7^) = 2j(sf + s 2 ). In this situation, an effective strong coupling 
is achieved for g c s > (7^, k). In the following experiments, this condition is fulfilled for weak 
external pumping at a fixed detuning of the drive laser (Ai/2n ~ 350 MHz). Note that this laser 
also provides weak Doppler cooling to reduce the motion of the ion in the trappy. 

We performed numerical simulations using a master equation approach to predict measured 
observables, i.e. the intra-cavity mean photon number, n ss , the Mandel Q-parameteP and the 
ratio R between coherent and incoherent emission rates on the lasing transition for continuous 
application of all lasers. The Q-parameter, Q = n ss (g^ 2 \Q) — 1), quantifies intensity fluctuations 
in the quantum (Q < 0) and classical (Q > 0) domains. The i?-parameter is defined as R = k ■ 
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n ss / (72 -ftp), where np is the steady-state population of the Pi/ 2 level, which controls the amplitude 
of incoherent photon scattering. Experimentally, R is accessed by monitoring the fluorescence rate 
on the S 1/2 P 1/2 transition which leads to np. 

In the following, we show that the ion-cavity device acts as a thresholdless laser when op- 
erated at the onset of the strong coupling regime, i.e. for g cS ~ 7^. Raman resonance be- 
tween selected Si/ 2 and D 3 / 2 states is established at low recycling rates by adjusting the drive 
laser detuning 22 . In this situation, the strength of external pumping is given by the recycling laser 
parameter, s 2 , which is held at a fixed detuning A 2 / (2n) = —20 MHz. 

Fig. 2 shows the measured steady-state photon number n ss and the Q- and ^-parameters as 
a function of external pumping together with simulation results (see supplementary information) 
for continuous laser excitation of the trapped ion at the boundary of the strong coupling regime 
(<? e ff ~ 7eff )• The following three regions are identified: At low recycling parameters (region 
(I) of Fig. 2), we observe the signature of a thresholdless single-atom laser. It is marked by an 
increase of the intra-cavity mean photon number, n ss , while intensity fluctuations, as indicated 
by the Q-parameter, are reduced. Lasing occurs in the quantum domain (Q = —2.3(4)% for 
n 2 /(27r) ~ 4.6 MHz) since single photons are stimulated into the cavity mode via the vacuum 
field while spontaneous emission is reduced (R m 0.7 . . . 0.8). In region (II) of Fig. 2, the laser 
starts losing its quantum properties as it approaches the classical domain (Q > 0), while the intra- 
cavity mean photon number reaches a maximum. For even stronger pumping, the recycling laser 
introduces phase noise leading to an increase of the Q-parameter (region (III) of Fig. 2). At the 
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same time it imposes a Stark shift of the D 3 / 2 levels, effectively reducing the atom-cavity coupling 
and therefore the mean photon number. This is a manifestation of the theoretically predicted self- 
quenching of single-atom lasers^^l 

In Fig. 3a-c, we have investigated the dramatic change in the statistics of the emitted light 
observed in regions (I)-(III) of Fig. 2 in more detail for a slightly different (albeit compatible) set of 
experimental parameters. Fig. 3a shows the measured and calculated normalized correlation func- 
tions g( 2 \r) for a small total cycling rate r tot /27r ~ 100 kHz. The slow recycling process intro- 
duces a time gap between successive photon emissions and sub-Poissonian light (Q ~ —1.2(1) %) 
is emittecP-^l Different timescales are apparent in Fig. 3a: The decrease of the correlation func- 
tion for | r | < 1.5 /xs reveals the finite storage time for cavity photons. For longer times r, the 
-function approaches unity on a timescale given by the total cycling rate r to t- Finally, for long 
time intervals r, a modulation of the correlation function is observed which can be attributed to 
the beat note between the ion's two quasi-degenerate radial modes of motion in the trappy. As 
we slightly increase the recycling parameter, the photon flux remains almost constant, but the sta- 
tistical properties of the emitted photons change dramatically: The correlation function becomes 
almost flat [Q = 0.0(1) % in Fig. 3b], akin to a Poissonian source of light. At even higher recycling 
parameters, the output statistics resembles that of a thermal source [Q = +1.0(1) % in Fig. 3c]. 
In general, the measured coincident photon correlations, (0), are governed by the population of 
the cavity Fock states with two photons, p 2 , so that g^(0) ~ 2p 2 /p\, where p\ is the population 
of states with only one photon. The population p 2 is established through an interplay between the 
rate r tot at which photons are added to the cavity and the rate 2k at which the photon state leaves 
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the cavity (see supplementary Table 2). The tunability of this population p 2 and thus the statistical 
properties of the emitted light is a feature of our system, which can be adjusted to operate in the 
strong or weak coupling regimes, as in the following where we demonstrate a single-atom laser 
exhibiting threshold behavior. 

In an experiment with different parameters, we studied the single-ion device for an even 
larger number of photons in the cavity. This was achieved by increasing the drive laser intensity to 
fli= 2ir ■ 130 MHz. Experimental results are displayed in Fig. 4 where we identify the following 
regimes: For low recycling rates, i.e. in region (I), the photon number and Q-parameter both rise 
indicating classical non-lasing behaviour. Increasing the external pumping further (region (II)), n ss 
still grows while the Q-parameter has passed through a maximum and decreases. Indeed, we note 
a 17 % increase of the intra-cavity photon number together with a reduction of 33 % in intensity 
fluctuations. At the same time, the rate of coherent vs. incoherent emission on the lasing transition 
indicated by the ratio R is on the order of unity. These are the signatures for the onset of a laser 
thresholdP-^2 As for a conventional laser, the threshold is reached when the population of photons 
in the cavity exceeds a critical value for which photon stimulated emission compensates incoherent 
processes. Most notably, this critical value is less than one for our one-atom device. Furthermore, 
we note that the variation of n ss at threshold is less pronounced than for classical lasers. This is a 
consequence of the low photon number in the lasing mode, as theoretically predicted HKKlGmil 
caused by a saturation of the rate of external pumping (r tot ) in region (II). Additionally, we observe 
in region (III) quenching of the laser output. Unlike the experiments presented in Fig. 2, this is 
not due to a level Stark shift since the Raman detuning is chosen such that resonant coupling is 
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established at threshold. Rather, numerical simulations indicate that it is due to the formation of 
dark states by coherent coupling via £l 2 , which disables the recycling process. More precisely, the 
induced dressed states, generally composed of two D 3 / 2 magnetic sub states and one Px/ 2 sub level, 
lose their P 1( / 2 contribution in region (III), where fi 2 is large. This renders the external pumping 
inefficient hence yielding self-quenching which is a striking difference compared to conventional 
lasers. In fact, it constitutes an intrinsic limitation of single-atom lasers. In Fig. 4, we also 
present theoretical simulations which agree qualitatively with the experimental data. We attribute 
the discrepancy to the non-zero temperature of the ion, which modifies the excitation spectrum for 
the parameter range of the experiments presented in Fig. 4 and is not included in the simulations 
(see supplementary information). 

In conclusion, we have demonstrated that a laser threshold can be achieved at the single-atom 
level, although much less pronounced compared to classical lasers^^ 12 * 16 * 17 ! as a consequence of 
the low photon number in the lasing mode. Another striking difference between one-atom and 
classical lasers is the ability to tune the photon emission statistics from quantum to classical. Such 
a light source offers perspectives to study novel types of spectroscopy, based on narrow band 
optical excitation with quantum/classical statistical distributions^. Our trapped ion laser is ideally 
suited to further investigate the transition between quantum and classical lasers through controlled 
addition of more and more ions interacting with the light field. 

Supplementary Information is linked to the online version of the paper at www.nature.com/nphys. 
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Figure 1. Experimental setup and excitation scheme, a, Schematic experimental setup: The ion 
(blue ball) is surrounded by a high-finesse cavity and positioned in an anti-node of the standing 
wave cavity field. The single-ion device is excited to a meta-stable state via a Raman transition 
(blue drive laser beam and cavity). A red laser provides repumping to close the excitation cycle. 
Photons leaving the cavity are coupled into a multi-mode optical fiber, split by a non-polarizing 
beam splitter and detected by two avalanche photodiodes (APD) in a Hanbury-Brown&Twiss 
setup. Photon arrival times are recorded with ps 300 ps temporal resolution, b, Relevant level- 
scheme of the 40 Ca + ion. A quantization axis is defined by applying a weak magnetic field, 
B = 2.8 G, perpendicular to the cavity axis. The drive laser (f^) is polarized parallel to the B 
field: It excites transitions conserving angular momentum (Am = 0). The cavity (g) supports two 
linear polarizations, perpendicular and parallel to the B field. The recycling laser (fi 2 ) is polarized 
perpendicular to the B field. 
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Figure 2. Thresholdless single-atom laser. Statistical behavior of the single-ion device at the 
boundary of the strong-coupling regime. Shown are the mean intra-cavity photon number, n ss , 
the Mandel Q parameter and the ratio R between coherent and incoherent scattering. Region (I) 
shows thresholdless lasing of the single-ion device, marked by decreasing intensity fluctuations 
(~ Q) and increasing photon flux (~ n ss ) as a function of external pumping, while R is on the 
order of one. In regions (II) and (III) the laser turns off (self-quenching). Drive laser parameters 
are (Q lt A t )= 2tt • (88, -350) MHz leading to W(2tt) = 165 kHz ~ ^/(2n) = 170 kHz. The 
recycling rate is varied by the Rabi frequency, Q 2 , for a constant detuning A 2 /2n w —20 MHz. 
Experimental data are obtained with the same trapped ion continuously probed over 12 hours. 
Solid lines show the results of numerical simulations using independently calibrated experimental 
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parameters. Error bars are derived from calibration uncertainties of the corresponding quantities 
(see supplementary information). 
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Figure 3. Tunable photon statistics. Normalized second order photon correlations, g^ 2 \r) for 
different photon generation rates. Increasing the photon generation rate leads to a transition from 
sub-Poissonian (a) via almost Poissonian (b) to super-Poissonian, bunched (c) light emitted by the 
cavity. The peak observed around r = /xs is a feature of the intermediate coupling regime and 
stems from the accumulation of photons inside the cavity. The drive laser excitation is such that 
(Q u Ai) « 2?r(95, -400) MHz while the recycling laser is set at (fi 2 , A 2 ) « 2?r(7, -20) MHz 
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in a, at (Q 2 , A 2 ) « 2vr(12, -20) MHz in b, and at (Q 2 , A 2 ) ~ 2tt(16, -20) MHz in c. In every 
measurement the cavity detuning is A c /2ir —400 MHz, the time resolution is set at 500 ns and 
accidental correlations have been subtracted (see supplementary information). Solid lines show 
the results of numerical simulations using independently calibrated experimental parameters. 
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Figure 4. Comparison between experimental observations (a) and theoretical predictions (b) for 
the statistical response of the single-ion device operating as a laser with threshold behaviour. 
Shown are the mean intra-cavity photon number, n ss , the Mandel Q parameter and the ratio R 
between coherent and incoherent scattering. Intensity (~ n ss ) and intensity fluctuations (~ Q) 
rise as a function of external pumping in region (I), indicating classical behaviour. Signatures of 
a laser at threshold are observed in region (II) where Q has passed through a maximum while n ss 
is still rising and R is on the order of one. Self-quenching turns the laser off in region (III). Drive 
laser parameters are (f2i, Ai)= 2n- (130, -350) MHz, resulting in a g c s ~ 240 kHz. The recycling 
laser has a fixed detuning A 2 /2n w —20 MHz. Experimental data are obtained with the same 
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trapped ion continuously probed over 12 hours. The dashed lines are a guide to the eye. Error 
bars are derived from calibration uncertainties of the corresponding quantities (see supplementary 
information). 
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1 Calibration procedures 

(N 
O 
O 

Cavity: The cavity decay rate was determined using a ring-down measuremenP with an error of 
less than 5 %. We measured the frequency spacing between the TEM 00 and TEM i modePto 
obtain the exact cavity length. This allows us to calculate the cavity mode volume, and thus the 
ideal ion-cavity coupling rate, g/2n = 1.61(2) MHz 3 . However, we estimate that the ion-cavity 
coupling reduces to g/2n = 1.3(1) MHz due to residual motion of the ion with respect to the trarP. 
In fact, our experiments are carried out at a large detuning of the drive laser, Ai/27r w 350 MHz, 
for which laser cooling has a reduced efficiency. 
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The length of the cavity is stabilized by servoing a piezo-electric transducer (PZT) attached to 
one of the cavity mirrors. The stabilization signal is derived using a transfer- lock scheme involving 
a frequency-stabilized laser at a wavelength of 785 nm which is resonant with one of the cavity 
modes. We estimate the stability of the cavity resonance frequency to be on the order of a few 
hundred kHz on a timescale of milliseconds. The instantaneous linewidth is expected to be close 
to the natural linewidth of the cavity. The relative position between the ion and the cavity can be 
adjusted by biasing the PZT attached to the second mirror of the cavity. 

We observe the cavity excitation spectrum by applying a constant (Rabi frequency and de- 
tuning) drive laser while sweeping the cavity detuning across the Raman resonance by tuning the 
785 nm laser frequency. Choosing a specific Raman resonance and fixed laser intensities and de- 
tunings, the cavity standing-wave can be mapped by recording the intra-cavity photon number via 
photon leakage through one of the cavity mirrors^ as a function of the bias voltage on the second 
mirror of the cavity. From the visibility of the standing wave, the localization of the ion relative to 
the cavity can be inferred to be 100 nrrP^. Derealization is attributed to two effects with approxi- 
mately equal contributions: (i) residual motion of the ion in the trap and (ii) mechanical vibration 
of cavity with respect to the trap. 

Magnetic field and trap frequencies: A bias magnetic field is induced by a pair of near-Helmholtz 
coils. The magnitude of the magnetic field was first measured by resolving the frequency of several 
Zeeman lines on the quadrupole transition Si/ 2 -H- D 5 / 2 at 729 nm 7 . This procedure was repeated 
for different currents, establishing a calibration curve with less than 5 % error. Additionally, the 
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trap frequencies were determined by identifying the motional sidebands on the Zeeman lineP. 

Laser parameters: The relevant laser parameters are calibrated in the absence of ion-cavity inter- 
actions. We record a collection of excitation spectra of resonance fluorescence to determine and/or 
set the excitation laser parameters, (0^2, Ai )2 ) and the laser linewidths 5i )2 of the drive (index 
1) and recycle (index 2) lasers. For each spectrum, the light emitted by the ion on the Su 2 -Pi/2 
transition is monitored while the recycling laser detuning, A 2 , is scanned. For the calibrations we 
typically record between five and ten spectra for fixed Cli and Ai while f2 2 is varied between each 
spectrum. Similar spectra are recorded for different settings of fii and Ai. Fitting all acquired 
spectra to a numerical model containing all 8 atomic states provides us with a calibration of all rel- 
evant experimental laser parameters with uncertainties as shown in Supplementary Table 1. Please 
note that all laser intensities are actively stabilized and the stability of our experimental conditions 
is tested by periodically (approximately every 30 min) probing the ion. The position of the ion in 
the cavity standing wave as well as possible laser frequency drifts are then continuously optimized 
during our experiments. 
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parameter 


typical value [MHz] 


relative error 


Ax 


350 


2.5 x 10~ 3 




... 250 


5% 


*i 


0.03 


12% 


A 2 


-20 . . . 20 


10% 




... 80 


5 % 




0.2 


70% 




350 


2.5 x 10~ 3 



Supplementary Table 1: Estimated uncertainties for the calibration of relevant experimental laser 
parameters. Note that the linewidth 5 2 of the recycling laser has negligible influence on the exper- 
imental results. 

Detection: We use two silicon avalanche photodiodes (PerkinElmer SPCM-AQR-15) with a mea- 
sured quantum efficiency of « 42% at 866 nm and a specified dark count-rate of < 50 counts/s in a 
Hanbury-Brown&Twiss setup to detect photons emitted by the cavity (Figure 1). The arrival times 
of detected photons are then time-tagged (PicoQuant PicoHarp 300) with a sub-nanosecond time 
resolution. At steady-state, including all detection losses of our experimental setup, an intra-cavity 
mean photon number of one corresponds to a photo-electronic count-rate of 32(4) kHz. 

The photons spontaneously emitted on the S1/2-P1/2 transition are collected by an objective 
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and are detected using a photomultiplier tube. The photon detection efficiency of this channel 
was determined from the fitting of a large series of excitation spectra of resonance fluorescence: 
7/397 = 1.8(1)%. Hence, the steady-state population of the Pi/ 2 levels is inferred, as well as the 
value of the i?-parameter. 

2 Evaluation of the - correlation function 

The photons emitted by the cavity are time-tagged with high efficiency by two fast photodiodes 
as described above. The detection events are then grouped into 500 ns bins and the g ^ -function, 
i.e. second-order time correlations among detection events, is calculated. However, the accidental 
correlations must be subtracted. These arise when background photons (stray light photons or 
APD dark counts) are correlated with themselves or output photons of the cavity. To quantify the 
amount of accidental correlations we use the signal to noise ratio (S/N) from the cavity excitation 
spectrum measured together with the correlation function. Considering that the noise follows a 
Poissonian time-distribution, the normalized correlation function can be calculated from the raw 
second-order correlation function of the detection events (#s+n( t ))' an ^ reads 8 



3 Theoretical model 

General description: Our theoretical predictions are based on a master equation approach for 
the total density matrix of the coupled ion-cavity system, p. The latter is expanded including eight 
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(1) 



relevant electronic levels for the ion part and two degenerate cavity modes with crossed polarization 
(h and v), in a number basis with at most four Fock states for each mode. The laser excitation 
of the ion is treated semi-classically whereas the ion-cavity interaction is written in the Jaynes- 
Cummings form in the rotating wave approximation!^. All decays, losses and laser linewidths have 
corresponding Liouvillians entering the master equation expressed in the Lindblad form. The time 
evolution of the density matrix, p{t), and the steady-state p ss are hence obtained by numerically 
solving a master equation for our experimental situation. 

Steady-state photon number: The steady-state intra-cavity photon number is directly obtained 
from p ss . It reads n ss =(n) ss , where n = (a\ah + ala v ), and a\ and a\ are the creation operators 
for cavity photons with horizontal (along the atomic quantization axis) and vertical polarization, 
respectively. Here we have defined the expectation value of the photon number, n, with respect 
to the steady-state density matrix, i.e. (n) ss =Tr[np ss ], Tr[...] denoting the evaluation of the trace 
along the ion-cavity subspace. 

Second-order correlation function (r) and Mandel Q -parameter: The second-order corre- 
lation function, g( 2 > (t, t + r), gives the probability to detect a photon at a time (t + r) conditioned 
on a previous photon detection at time t. It is usually referred to as intensity correlation function 
and reads 

m, s (■ I(t)I(t + T) :) 

where the intensity-operator, I(t), is proportional to the intra-cavity photon number, I(t) oc n(t). 
Furthermore, note that operators inside colons must be time and normally ordered in Supplemen- 
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tary Equation (O. At steady-state, i.e. for t — > oo, the normalized correlation-function is a function 
of r only and it follows 

9 < 2 » (r) = Wl>>£. (3) 

In Supplementary Equation ©, the expectation value (...) is evaluated using the conditional den- 
sity matrix p = {a-hPss^h + a vP ss o.l)/n ss . It shows that the (/^-function resolves the transient 
evolution of the intra-cavity photon number, upon annihilation of a photon at steady state. It di- 
rectly relates to intensity fluctuations via the Mandel Q-parameter which reads Q=n ss (g( 2 \0) — 
l)=(n) - n ss . 

Finally, let us mention that our theoretical model allows us to evaluate any observable of the 
single-ion device. For instance, we confirm in the following that the measured correlations for low 
recycling rates g^ 2 \0) ~ 2p 2 /Pi (Figure 3) by extracting the populations p\ and p 2 for the one- 
and two-photon Fock state, respectively. Supplementary Table 2 compares our simulations with 
the experimentally obtained 5^(0). 
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Supplementary Table 2: Simulated one- and two-photon Fock state populations. This leads to 
g^(Q) ~ 2p 2 /Pi which is compared to the experimentally observed g^(0). 
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Influence of the motion of the ion: For the experiments reported in Figure 2, we reach a quan- 
titative agreement between our observations and the results of our theoretical simulations solely 
considering independently measured experimental parameters. By contrast, for stronger excitation 
of the drive laser, as in Figure 4, we achieve only qualitative agreement. 

In Figure 4, theoretical simulations are compared to experimental findings when considering 
calibrated experimental parameters. Qualitatively, these exhibit the same statistical variations as 
a function of external pumping, but the measured intra-cavity mean photon number differs from 
the predictions. Consequently, the amplitudes of Q and R differ as well. We mainly attribute 
the observed discrepancy to the motional degrees of freedom of the ion which are not considered 
in our model. In the pumping process, the drive laser for the Raman transition provides only 
weak cooling, since it is far detuned from the S1/2- P1/2 resonance. This results in a non-negligible 
motional excitation of the ion in the trap. As a consequence, the atom-cavity coupling, g, is reduced 
by up to 20% and motional sidebands involving the absorption or emission of phonons complicate 
the observed Raman spectrum 4 . 

In the experiments reported in Figure 4, the cavity mode is set to be on the slope of the Ra- 
man resonance at low recycling rates. Resonant coupling is then achieved in region (II) due to the 
Stark shift induced by the recycling laser. At the same time, the shape of the resonance and the 
resultant experimental parameters (n 8B , Q, R) are strongly distorted by the motional sidebandP. 
Thereby, quantitative agreement between simulation and experiment is lost. On the other hand, in 
the experiments presented in Figure 2, the single-ion device is operated close to the Raman reso- 
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nance throughout the investigated parameter range. Hence, the ion motional degrees of freedom 
play a minor role and theoretical predictions quantitatively reproduce experimental findings. 
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